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I.  INTRODUCTION 


Radio  comniunication  usage  continues  to  expand  because  the  military  and 
commercial  markets  seek  (1)  personalized  service  such  as  paging  and  cellular 
applications,  (2)  ever  greater  exchanges  of  text  and  data  (facsimile,  computer 
networking),  and  (3)  transfer  of  high  quality  signals  (television).  The  messages  are 
"traffic".  The  carriers  are  sine  waves. 

There  are  two  primary  ways  of  carrying  messages  using  sinusoids.  The  message 
can  cause  to  vary  (modulate)  the  amplitude  or  the  argument  (angle)  of  the  sinusoid. 
These  are  called  amplitude  modulation  (AM)  and  angle  modulation  (4-M) 
respectively. 

Amplitude  modulation  can  assume  various  forms  such  as  double  sideband  with 
carrier,  double  sideband  suppressed  carrier  and  single  sideband.  Angle  modulation 
can  be  either  frequency  modulation  (FM)  or  phase  modulation  (PM).  With  analog 
messages,  there  is  little  difference  between  FM  and  PM.  In  fact,  FM  hardware 
would  be  used  should  PM  be  chosen.  Consequently,  PM  with  analog  messages  is,  in 
practice,  reduced  to  FM.  [  Ref.  1:  pp.  198-199  ] 

Frequency  modulation  exhibits  several  interesting  characteristics  not  shared  by 
AM.  Two  such  characteristics  are  threshold  effect  and  capture  effect.  The  threshold 
effect  relates  to  the  quality  (signal-to-noise  ratio)  of  the  demodulator  output  as  a 
function  of  the  quality  of  the  received  signal  (demodulator  input).  Below 
"threshold"  quality  of  the  input,  the  output  quality  deteriorates  rapidly.  The 
capture  effect  relates  to  the  ability  of  the  demodulator  to  recover  the  message  of  the 
dominant  carrier  when  two  or  more  FM  carriers  are  present. 
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This  research  is  concerned  with  the  capture  effect  of  FM.  Computer  simulation 
is  used  to  verify  an  idea  concerning  the  reason  for  (origin  of)  the  capture  effect  of 
FM  demodulators. 

Background  information  on  FM  and  the  capture  effect  is  presented  in  chapter  II 
of  this  thesis.  Chapter  III  describes  the  computer  simulation  procedure  and 
program.  Results  of  the  simulation  are  contained  in  Chapter  IV.  Conclusions  and 
reconunendations  are  listed  in  chapter  V.  The  computer  program  is  listed  in  the 
Appendix.  A  list  of  references  is  provided. 
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11.  BACKGROUND 


A.  EQUATION  FOR  FM  CARRIER 

An  angle  modulated  sinusoid  is  expressed  as 

s{0  =  Bcos[<i»jO)]  (2.1) 

where  the  constant  B  is  the  carrier  peak  amplitude.  A  complete  oscillation  cycle 
occurs  whenever  the  argument  ^{t)  changes  by  2%  radians.  The  instantaneous 
frequency  of  the  angle  modulated  wave  s(t)  is  defined  by 


~Si 


(2.2) 


Frequency  modulation  is  that  form  of  angle  modulation  for  which  the 
instantaneous  frequency  f-{t)  varies  in  accordance  with  the  baseband  signal  or 
message  voltage  m{t).  Typically,  a  linear  relationship  is  chosen  so  that 


fft)  =  fc+ (2.3) 
where  is  the  frequency  of  the  unmodulated  carrier  and  the  constant  kj  represents 
the  frequency  sensitivity  (in  Hertz  per  volt)  of  the  modulator.  From  Eqs.  (2.2)  and 
(2.3),  the  angle  (pit)  is 

(p^{t)  =  2tf^t  +  2wkjf^m{T)dT  (2.4) 

where,  for  convenience,  the  angle  of  the  unmodulated  carrier  wave  is  assumed  to  be 
0  at  t  =  0.  The  frequency  modulated  wave  is,  therefore, 

s(t)  =  B  cos[25r/ 1  +  2Ttkff^m{T)dT\  (2.5) 
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B.  CAPTURE  EFFECT 

Consider  the  case  of  two  FM  carriers,  si(<)  and  S2{t),  present  at  the  demodulator 
input  and  in  the  same  frequency  band.  It  is  well  known  in  practice  that  the  output 
of  the  demodulator  is  the  message  of  the  dominant  carrier.  That  is,  if 

Si(t)  =  Bi  cos[27rf  t  +  27rkf  f^mi{T)dr] 
c  JJq 

S2{t)  =  B2  cos[2ff(/+  i)t  +  2tkf  f^m2{T)dT]  (2.6) 

C  JJq 

then  the  output  of  the  demodulator  is  proportional  to 

mi(t),  when  B\>B2 
Tn2it),  when  Bi<B2 

The  carrier  frequency  offset  t  is  any  value  such  that  both  si(  t)  and  S2{  t)  are  in 
the  same  operating  frequency  band  of  the  receiver/demodulator.  This  inherent 
ability  of  an  FM  demodulator  to  suppress  weaker  FM  carriers  is  called  capture 
effect. 

The  FM  demodulator  can  be  a  phase  locked  loop,  a  slope  detector  which 
includes  Foster-Seeley  and  ratio  detectors,  or  a  pulse  counting  discriminator. 

C.  ANALYSIS 

As  shown  by  Eq.  (2.3),  the  instantaneous  frequency  fit)  of  a  frequency 
modulated  carrier  is  proportional  to  the  message  m{t)  being  carried.  Consequently, 
it  is  claimed  in  the  literature  that  the  output  voltage  of  a  frequency  demodulator  is 
proportional  to  fU)  —  L  =  A /•(<)  [  Ref.  2:  pp.  169—172  ].  This  is  indeed  the  case  in 
practice  when  the  demodulator  input  is  a  single  frequency  modulated  carrier. 
However,  it  is  not  the  case,  because  of  the  capture  effect,  when  the  demodulator 
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input  is  multiple  frequency  modulated  carriers. 

The  instantaneous  frequency  of  the  sum  of  two  frequency  modulated  carriers 
can  be  derived  as  follows.  L„t 


v(t)  =  Bi  cos^i(t)  +  Bi  cos(^(t) 
Using  trigonometric  identities,  Eq.  (2.7)  can  be  written  as 

u(i)  =  Cco8(t>z{t) 


where, 


C\i)  =  B,2  +  B22  +  2BiB2CO^\<i>i{i)-(h{i)] 

Bis\n<^i{t)  +  B2s\n<p2{t) 

tan(^3(<)  = - 

BiCOs4>tit)  +  B2COS  (hit) 

From  Eq.  (2.2),  the  instantaneous  frequency  f^t)  is 


(2.7) 


(2.8) 


(2.9) 

(2.10) 


where, 


1  d 
U 


tan"' 


’5isin<^if< 

+B2Sin(/>2(U' 

filCOS0i(t 

+  B2COS(^(t) 

_  (Bi</>i'cos(/>i+B2<^'cos<^2 )  C+(^i(;^i'sin(?^i+B2(/>2'sin(fe)Z? 

27r(C2^./)2) 

_  B\^d)\  -\-B2^(i>2  +  [  ^1^2COS(^1-02)1(  <i>\—<h') 

2  7r[fii2+fi2^'l'25i  B2C0S(  ^J“^2  )] 


(2.11) 


C=  B\  zos(f>\{t)+B2  cos<^(<) 
D=  Bi  sin(?^(<)+^2  sin02(t) 
and  0i',  <p2'  represent  and  respectively. 


5 


As  before,  for  an  FM  carrier 


4>-{t)  =  2;r/^  +  2irk^m^{t)  i  =  1,  2  (2.12) 

where  the  carrier  frequencies  for  both  messages  are  assumed  to  be  the  same  and 
m  / 1)  =  i  th  message 

I 

=  frequency  sensitivity  of  the  frequency  modulator  in  Hertz/ volt 

Equations  (2.8)  through  (2.12)  can  be  reduced  and  interpreted  for  particular 
messages.  For  example,  let  the  two  messages  be  tones  where 

m-(<)  =  D.  cos2x^<t)  i  =  1,  2  (2.13) 

For  simplicity,  let  A  =  A  =  1-  Then, 

=  2;r/J  +  2tkf  f^cos{2nfT)dT 
kj 

=  2nfJ  +  -sl-^2irf.t 

h 

<f>'  (t)  =  2x/  +  27r^rCOS2x/.i  j  =  1,  2  (2.14) 

I  C  j  X 

Using  these  values  in  Eq.(2.11)  gives 


where 


UBiHB2^)+kf'X+BiB2Y^Z 

f.{t)  =  S. - 1 - 


(2.15) 


X  =  cos27r^t  +  cos27rj^t 
Y  =  2f^+k j{cos2nJit  +  cos2ir^t) 

Z  =  cos[2;rA:y(  sin2?r/2<  +  sin2?r/i^ )] 


I 
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Equation  (2.15)  is  plotted  in  Figs.  2.1.a  through  2.2.c  along  n-ith  mi{t)  and 


m2{t)  where  /i  =  1  IIz  and  /2  =  1-3  Hz.  It  can  oe  seen  that  there  are  rriany 
impulselike  frequency  deviations  from  the  message  waveforms.  The  variable  in  Fig. 
2.1  is  /?  (  modulation  index  for  FM).  The  bandwidth  of  the  modulated  carrier  is 
related  to  /?.  In  general,  bandwith  increases  linearly  with  0  [  Ref.  3:  pp.  194—197]. 

There  are  two  major  factors  that  affect  the  graph  of  instantaneous  frequency. 
First,  the  bandwidth  of  the  modulated  carrier  affects  the  spacing  in  time  between 
adjacent  impulse— like  frequency  deviations.  Figures  2.1. a  through  2.1.C  show  that  as 
the  bandwidth  decreases,  the  spacing  increases.  Second,  the  ditference  in  magnitude 
of  the  two  modulated  carriers  affects  the  extent  of  the  deviation.  Figures  2. 2. a 
through  2.2.C  show  that  as  the  magnitude  difference  increases,  the  deviation  extent 
decreases. 

However,  these  results  do  not  show  specific  evidence  of  capture.  Rather,  they 
suggest  that  the  instantaneous  frequency  needs  to  be  averaged  over  a  short  time 
interval  to  reveal  the  capture  effect.  The  same  thing  is  true  for  the  case  of  close  (but 
different)  carrier  frequencies,  as  shown  in  Chapter  III. 

The  idea  that  the  capture  effect  is  caused  by  the  lowpass  filtering  operation  of 
frequency  demodulation  led  to  this  research.  Since  lowpass  filtering  can  be  thought 
of  as  averaging  (or  smoothing),  it  was  decided  to  investigate  the  capture  effect  by 
averaging  the  instantaneous  frequency  using  a  computer. 

So,  in  hardware,  the  averaging  is  done  by  the  lowpass  filter  (LPF).  A  LPF  is 
part  of  every  frequency  demoaulator  (fiequency—to— voltage  or  voltage— to-frequency 
converter).  Consequently,  the  output  of  a  frequency  demodulator  is  a  lowpass 
filtered  version  of  instantaneous  frequency  deviation  A/^t). 
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Again,  lowpass  filtering  can  be  considered  as  a  smoothing  operation  on  the 
applied  voltage.  In  this  research,  we  used  computer  simulation  to  investigate  the 
effect  of  smoothing  of  ^f^t)  to  reveal  the  presence  of  capture.  Smoothing  is  created 
by  averaging  time  intervals  between  adjacent  zero  crossings  of  the  FM  carrier.  Since 
this  interval  ( Ti)  represents  the  time  required  for  t  radians  excursion  of  the 
argument  of  the  cosine  carrier,  and  since  t  radians  is  taken  as  a  half  cycle,  then  two 
such  intervals  (Ti  +  T^)  become  a  full  cycle  and  number  of  cycles  per 

second  or  the  instantaneous  frequency  for  7i  +  T-i  seconds.  Averaging  N  such 
groupings  provides  a  smoothed  value  of  or  short-term  average  of  the  instantaneous 
frequency  over  that  expanded  time  interval. 

A  computer  program  was  used  to  approximate  T\  and  for  every  interval 
between  zeros  of  v{t).  These  values  are  then  averaged  and  the  results  are  plotted. 
The  program  is  explained  in  Chapter  III. 
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Frequancy 


Fr*qu«neu  Frequency 


Fraquancg 


Fig.  2.2  Instantaneous  Frequency  for  Tone  Messages 
where  0  =  10,  =  200  Hz,  /i  =  1  Hz,  =  1.3  Hz 
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m.  SIMULATION  PROCEDURE 

The  FM  demodulator,  including  LPF,  is  simulated  by  a  computer  program  so 
that  critical  parameters,  such  as  signal  amplitudes,  bandwidth  of  the  FM  system, 
message  forms,  and  interval  of  averaging,  can  be  easily  changed.  Since  the  sampled 
version  of  the  demodulator  input  forms  a  sequence  of  numbers,  it  is  convenient  to 
consider  it  as  a  row  (or  column)  vector.  Once  that  is  done,  the  rest  is  simple 
manipulation  of  the  row  (or  column)  vector. 

A.  WAVEFORM  GENERATION 

The  desired  waveform  v{t)  can  be  generated  using  a  simple  routine 
(WAVEGEN).  In  this  routine,  the  carrier  amplitudes,  values  of  kp  carrier  frequ¬ 
encies  and  message  frequencies  are  input  parameters  which  can  be  changed.  This 
resulting  wave  is  sampled  at  a  rate  which  will  provide  the  desired  accuracy  in 
measuring  the  time  between  adjacent  zeros  of  v{t).  In  this  work,  a  carrier  frequency 
of  200  Hz  and  sampling  rate  of  12  kHz  are  used. 

The  sampled  values  contain  no  information  about  frequency.  Rather,  the 
number  of  samples  between  zero  crossings  are  related  to  frequency  because  the 
number  of  samples  between  zero  crossings  itself  represents  a  half  cycle  oscillation  of 
argument  of  u(<)  and  its  reciprocal  can  be  interpreted  as  twice  the  frequency  of  that 
period.  Therefore,  the  routine  notes  the  polarity  of  each  sample  so  that  the  number 
of  samples  between  zero  crossings  (polarity  changes)  can  be  easily  counted.  This 
result  is  then  used  as  the  input  for  the  next  routine. 
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n.  COUNTING 

The  number  of  samples  between  zero  crossings  (polarity  changes)  multiplied  by 
the  fixed  sampling  interval  (A)  is  considered  to  be  taken  as  half  the  instantaneous 
period  of  that  interval.  So,  a  sequence  of  integers  which  consists  of  the  number  of 
samples  between  zero  crossings  need  only  be  generated.  The  i  th  entry  for  this 
sequence  is  denoted  by  7'^- ,  j  =  1,  2,  3,  •  •  •  as  shown  in  Fig.  3.1.  The  totality  of 
entries  forms  a  sequence  of  integers. 


Fig.  3. 1  Formation  of  sequence 

Furthermore,  T-  itself  is  a  time  index.  Since  there  are  » -  1  entries  before  the 

j-1  . 

current  interval  T- ,  is  the  period  for  the  time  +  •  (The  middle  of 

the  interval  is  selected  for  the  time  index  of  interval  .  For  example,  the  time 
index  of  T3  is  7i  +  72  +  ^Ts  ) 

Since  the  first  and  last  interval  may  contain  a  number  of  samples  that  does  not 
correspond  to  a  complete  interval  between  adjacent  zeros,  these  two  entries  are  not 
used  in  the  averaging  routine. 
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C.  AVERAGING 


Given  the  sequence  of  integers,  any  average  over  N  intervals  is  formed  as 

k+N 
E  T, 
k  * 

T,  = -  k  =  any  positive  integer  (3.1) 

«  AT 

T-  =  I  th  interval 

t 

N  =  number  of  intervals  averaged 


Since  frequency  is  taken  as  the  reciprocal  of  the  period,  the  averaged  frequency, 
denoted  by  ^ ,  is  given  by 


1 

2T^  a 


(3.2) 


To  obtain  more  realistic  (smoother)  results,  a  moving  average  is  used  rather 
than  fixed  time  interval  averaging.  To  illustrate  this,  assume  a  sequence  of  intervals 
as  follows 


Ti  T2  n  Ti  n . 

Then,  if  the  average  is  done  over  3  intervals,  the  averaged  version  of  the  sequence 
will  have  the  following  entries 


T1+T2+T3 
- 3 - 

T2+r3+T4 

3 - 


'^k  '^k+1 


+  T 


Ar+2 


(3.3) 
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Again  from  the  relationship  between  period  and  frequency  (Eq.(3.2)),  averaged 

version  of  frequency  will  form  the  sequence 

_ 1  _ 1  1 

2T  X  2T2X  2T^.X 

In  this  case,  the  time  indices  for  the  averaged  version  of  intervals  are 

TH-Ta+Tj  ^ 

Ti+— —  =  time  index  for  Ta 

Ti+  =  time  index  for  ^3,  and  so  on. 

Generally,  time  indices  for  the  moving  average  over  N  intervals  are  found  as 
i— 1  .  i+A-l 

(tinae  index  for  !r^. )  (3.4) 

Since  the  values  of  each  Ti  are  not  the  same,  the  averaged  frequency  values  when 
plotted  will  not  be  equally  spaced. 
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IV.  SIMULATION  RESULTS 

With  an  average  of  2  intervals  (N  =  2),  the  resulting  plot  is  very  similar  to  that 
of  the  analytical  one.  That  is,  when  almost  no  averaging  is  done,  the  resulting  plot 
should  yield  that  of  the  analytical  one.  This  is  another  way  of  verifying  that  the 
simulation  program  is  running  properly.  As  the  number  N  of  average  intervals 
increases,  the  impulse— like  frequency  deviation  becomes  smaller  and  the  overall 
frequency  deviation  approaches  that  of  the  stronger  signal. 

Figures  4.1.a  through  4.5.C  are  the  simulation  results.  The  results  are  presented 
as  averaged  instantaneous  frequency  versus  time.  In  all  cases,  the  message  mi(t) 
corresponds  to  that  of  the  dominant  carrier.  The  following  variables  are  used  to  see 
the  effect  of  system  parameters. 

ai/a2  =  ratio  of  carrier  amplitudes  (1  through  1.5) 

/?  =  modulation  index  (1  for  narrowband,  10  for  wideband,  50 
for  very  wideband  FM) 

J[  =  message  1  frequency  (1  Hz  fur  sinusoidal  message) 

^  =  message  2  frequency  (1.3  Hz  for  sinusoidal  message) 
fc  =  carrier  frequency  (180  Hz  through  200  Hz) 
iV=  number  of  intervals  averaged  (2,  20,  40,  60,  80) 

Various  messages  were  used  in  the  simulation.  Figures  4.1. a  through  4.1.C  show 
the  results  when  the  messages  are  linear  intrapulse  FM  (up  CHIRP  and  down 
CHIRP).  Figures  4.2.a  and  4.2. b  are  for  exponentially  decaying  and  growing 
sinusoidal  messages.  Sinusoidal  messages  are  used  extensively  to  show  the  results  for 
various  correlations  of  the  listed  variables. 


Even  though  there  may  be  other  factors  which  influence  the  capture  effect,  only 
the  following  four  factors  are  investigated  in  this  research. 

A.  EFFECT  OF  N  (AVERAGING  TIME) 

Figures  4. 3. a  through  4.3.d  shows  the  effect  of  N.  The  other  factors  (ai/fl2,  A 
fi,  X)  remain  fixed.  This  factor  is  essentially  the  cause  of  the  capture  effect.  As  N 
increases,  the  instantaneous  frequency  increasingly  resembles  the  message  of  the 
stronger  signal. 

B.  EFFECT  OF  SIGNAL  AMPLITUDE  RATIO 

Figures  4. 4. a  through  4.4. d  show  the  effect  of  signal  amplitude  difference.  While 
the  other  factors  remain  fixed,  ai/a2  varies  from  1.01  to  2.0.  The  effect  of  signal 
amplitude  difference  is  very  interesting.  Even  with  very  small  signal  amplitude 
difference,  the  system  still  locks  onto  the  stronger  signal.  For  numerical  comparison, 
the  0.8  dB  difference  case  (Fig.  4.4.b)  and  1.6  dB  case  (Fig.  4.4. c)  show  little 
difference.  In  the  case  of  very  small  amplitude  difference,  more  averaging  is  needed 
to  see  the  capture  effect.  Figure  4.4. a  indicates  capture  for  a  carrier  separation  of 
0.17  dB  (ai/c2  =  1.02). 

C.  EFFECT  OF  MODULATION  INDEX 

Modulation  indicies  are  changed  from  1  to  50  in  Figs.  4. 5. a  through  4.5.c. 
Capture  occurs  for  both  narrowband  FM  and  wideband  FM.  Since  the  narrower 
band  FM  shows  fewer  impulse— like  frequency  deviations,  each  of  longer  duration 
than  those  of  wider  band  FM,  narrow  band  FM  needs  to  be  averaged  longer  (see 
Fig.  2.1). 
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D.  EFFECT  OF  CARRIER  FREQUENCY  DIFFERENCE 

Figures  4. 6. a  and  4.6. b  show  the  effect  of  carrier  frequency  difference  on 
capture.  It  appears  that  capture  becomes  more  complete  as  the  carrier  frequencies 
separate.  Obviously,  the  two  carriers  must  lie  in  the  passband  of  the  receiver 
(intermediate  frequency  amplifier)  if  they  are  to  be  present  at  the  input  to  the 
demodulator. 
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Frequency  Frequency 


C.  N=iO 


d.  /V=80 

Fig.  4.3  Average  Iiislanlaneoiis  Frequency  Deviation  for  Tone  Messages 
where  01/02  =1.2,0=  10,  =  200  IIz,  /i  =  I  IIz,  =  1.3  IIz 


4 


2 


Fig.  4.4  Average  Inslantaneous  Frequency  Deviation  for  Tone  Messages 
where  /?  =  10,  /  =  2U0  IIz,  /i  =  1  Hz,  =  1.3  Hz,  A  =  60 


FrBqusncy 


c.  /?=50 

Fig.  4.5  Average  Instantaneou.s  Frequency  Deviation  for  Tone  Messages 
where  a{jai  =  1.2,  =  200  IIz,  /i  =  1  Hz,  ^  =  1.3  Hz,  =  60 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  simulation  results  establish  that  the  lowpass  filtering  portion  of  frequency 
demodulation  accounts  for  the  capture  effect  of  FM  receivers.  Without  filtering,  the 
instantaneous  frequency  deviation  of  the  sum  of  two  carriers  has  an  impulse— like 
history  which  does  not  consistently  favor  the  message  of  the  dominant  carrier.  With 
averaging  or  smoothing  (lowpass  filtering)  of  the  instantaneous  frequency  deviation, 
capture  occurs. 

The  simulation  results  indicate  capture  occurs  for  both  narrowband  and 
wideband  FM.  Capture  is  evident  when  the  carriers  are  separated  in  amplitude  by 
less  than  0.2  dB.  The  effect  of  averaging  time  is  revealed  by  the  simulation  output. 
In  general,  longer  averaging  times  produce  smoother  approximation  to  the  message 
of  the  dominant  carrier,  as  expected,  and  also  show  the  capture  effect  more 
convincingly. 

It  is  recommended  that  the  capture  effect  be  verified  using  an  operating 
experimental  system  in  which  system  parameters  can  be  controlled  and  accurately 
measured. 
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APPENDIX 


Program  List 

This  appendix  contains  the  program  used  to  calculate  f-  defined  in  the  thesis. 
The  program  actually  determines  T-and  as  given  by  Eq.  (3.1). 

The  input  equation  is  si(t)  +  S2(t)  as  defined  by  Eq.2.6.  This  equation  is 
sampled  12,000  times  in  the  one  second  simulation  interval.  The  interval  between 
samples  is,  then,  1/12,000  =  83.3  /isec  =  A.  The  polarity  of  each  sample  value  is 
noted.  The  number  of  consecutive  sample  values  of  like  polarity  in  interval  i  is  . 

The  program  stores  each  value  of  and  forms  T^,  from  Eq.  (3.1)  for  various 
assigned  values  of  N.  Then  is  calculated  according  to  Eq.  (3.2).  Finally,  ^is 
plotted  for  each  time  index  of  Eq.  (3.4). 

The  program  contents  are  following; 

MAIN  ;  main  routine  for  the  simulation  of  capture  effect 

CHIRP  GEN  :  subfunction  routine  for  generation  of  CHIRPs 
WAVEGEN  :  subfunction  routine  for  generation  of  sinusoidal  waves 
»  COUNT  ;  subfunction  routine  for  counting  and  averaging 

SAMPLE  :  subfunction  routine  for  display 
The  program  was  written  in  APL  and  run  on  IBM  PC/AT  with  3  Mb  of  memory. 
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A.  MAIN 


V  Y«-MAIN  A;A;  Y  ;TMrSMP;TSLICE 

[I]  n 

12}  n  This  is  main  routine  for  the  capture  effect 
13]  n  simulation.  For  the  clality  of  logical  flow, 
t4]  n  each  simulation  step  is  divided  into  their 
[5]  n  own  routines. 

t6]  n  The  output  will  be, 

[7]  n  T:  time  index  for  fl  and  f2 

[8]  o  TIME:  time  index  for  aceraged  frequency 

t9]  n  (not  equally  spaced) 

[10]  n  F1,F2:  individual  instantaneous 

[II]  0  frequency 

[12]  n  RES:  averaged  version  of  frequency 

[13]  n 
[14  ] 

[15]  Genera  ting  waves....' 

[16] 

[17]  n  Generating  arbitrally  sinusoidal  waves 

[18] 

[19]  TMPSMPt-WAVEGEN  A 

[20]  []♦-' Counting.  ...  ' 

[21] 

[22]  n  Count  the  number  samples  between 

[23]  0  zero  crossings  and  form  a  resulting  sequency 

[24] 

[25]  Y^-COUNT  TMPSMP 

[26]  U«-' Sampling.  ...  ' 

[27]  t 

[28]  n  Since  the  numbers  of  elements  of  sampled  time  and 

[29]  n  message  are  too  many  for  plotting  purpose,  it  is 

[30]  n  sampled  without  missing  any  important  information. 

[31] 

[32]  T^S AMPLE  T 

[33]  F1<-SAMPLE  Fl 

[34]  F2<-SAMPLE  F2 

[35]  G<-'Done.' 

V 
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B.  CHIRPGEN 


V  V^-CHIRPGEN  A1;DI0;A1;A2;PH1;PH2;V;TMP1;TMP 

[1] 

[2]  n  This  routine  generates  two  CHIRP  signals,  namely  up 

[3]  0  and  down  CHIRP  to  test  the  program. 

[4]  fi  Equations  for  this  routine  are, 

[5]  R  F1=100T+150  for  up  CHIRP 

[6]  R  F2=250-100T  for  down  CHIRP 

[73  R  where  OsTsl 

[83 

[9  3  DIO^-O 
[103  A2<-1 

[113  TSLICE«-r(l. 3*0. 53x10000 

[12  3  T<-(tTSLICE3^TSLICE 

[13  3  D<-50 

[14  3  TMP«-(02  3xD 

[153  TMP<-TMPxT 

[16  3  TMP«-TMPxT 

[173  TMPl<-(O2)x(200-D) 

[18  3  TMP1<-TMP1XT 
[19  3  PH1«-TMP+TMP1 
[203  TMP«-(O23x(200  +  D) 

[213  TMP«-TMPxT 
[22  3  TMP1«-(02)xD 

[23  3  TMP14-TMP1XT 
[24  3  TMPl*-TMPlxT 
[25  3  PH2*-TMP-TMP1 
[26  3  TMP+-20PH1 

[27  3  TMP1<-20PH2  » 

[28  3  TMP<-AlxTMP 
[29  3  TMP14-A2XTMP1 
[30  3  V«-TMP+TMP1 
[313  V«-xV 
[32  3  Fl<-2xDxT 
[33  3  F14-F1+200-D 

[343  F2«-200+D 

[35  3  TMP«-2XD 
[36  3  TMP<-TMPxT 
[37  3  F2«-F2-TMP 
V 
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C.  WAVE6EM 


V  V<-WAVEGEN  A1 ; DIO ;  A1 ;  A2  ;  V;TMP1  ;TMP ;  VI ;  V2  ;  KFl ;  KF2 

[ 1 ]  ; FCl ; FC2 ; FMl ; FM2 ; PHI ; PH2 ; PI ; P2 

[2]  n 

[3]  fl  Routine  to  generate  pa.nples  of  sinusoid  waveforms 

[4]  n 

[5]  DI04-0 
16] 

[7]  0  Second  wave  amplitude.  This  routine  accepts  the  first 

[8]  n  signal  amplitude  as  an  input  argument 

[9] 

[10]  A2«-l 

[11]  P2<-o2 

[12]  TSLICE<-r(1.5*0.5)xl0000 

[13]  T*-(V.TSLICE)^TSLICE 

[14] 

[15]  o  Carrier  frequency,  frequency  sensitivity  and  message 

[16]  0  frequency 

[17] 

[18]  FCl«-200  O  FC2<-200 

[19]  KF14-10  O  KF2<-10 

[20]  FM1<-1  O  FM2«'1.3 

[21]  PHl«-0  O  PH2«'0 

[22]  TMP<-P2xFMl 

[23]  TMP<-TMPXT 

[24]  TMP4-TMP+PH1 

[25]  TMP<-loTMP 

[26]  TMP«-TMP-(loPHl) 

[  27  ]  TMP«-TMPxKFl-f  FMl 

[28]  TMP1<-P2xFC1xT  ^ 

[29]  TMP^-TMP+TMPl 

[30]  V14-20TMP 

[31]  Vl<-VlxAl 

[32]  TMP4-P2XFM2 

[33]  TMP«-TMPxT 

[34]  TMP«-TMP+PH2 

[35]  TMP4-10TMP 

[36]  TMP4-TMP- C  10PH2  ) 

[37]  TMP«-TMPXKF2+FM2 

[38]  TMP1*-P2xFC2 

[39]  TMPl*-TMPlxT 

[40]  TMP«-TMP+TMP1 

[41]  V2«-2otMP 

[42]  V2<-A2xV2 

[43]  V«-x(Vl+V2) 

[44]  TMP+-P2XFM1 

[45]  TMP«-TMPxT 

[46]  TMP<-TMP+PH1 

[47]  TMP«-KF1XC0S(TMP) 

[48]  F1«-FC1+TMP 

[49]  TMP<-P2xFM2 

[50]  TMP«-TMPxT 

[51]  TMP*-TMP+PH2 

[52]  TMP«-KF2xCOS(TMP) 

[53]  F2+-FC2+TMP 
V 


D.  COUNT 


V  y«-COUNT  V; C ; COUNTER ;INIT; I; LIMIT ;D;Y;TMP 

[1]  ;EXTIME1;EXTIME2;INTV 

[2] 

[3]  fiRoutine  to  form  half  period  sequence  and  average  it 
C4] 

[5]  o 

[6]  n  Set  the  initial  matricies  and  values 

[73  0 

[  8  ]  C<-^ 

[9]  COUNTER<-0 

[10]  INIT+-VCDIO] 

[11]  I«-DIO 

[12]  LIMIT4-(pV)+niO-l 

[13]  n 

[14]  o  Routine  to  count  number  of  samples  between 

[15]  0  Zero-crossings 

[16]  fi 

[17]  LOOPl; 

[18]  -♦(V[I]=INIT)/L00P2 

[19]  C«-C,  COUNTER 

[20]  INIT«-V[I] 

[21]  COUNTERS- 1 

[22]  •♦L00P3  . 

[23]  LOOP2 ! 

[24]  C0UNTER4-C0UNTER+ 1 

[25]  L00P3: 

[26]  I*-I  +  l 

[27]  -♦(IsLIMITJ/LOOPl 

[28]  C<-C,  COUNTER  i 

[29]  0 

[30]  0  Drop  the  incomplete  cycle  at  the  beginning  and  the  end 

[31]  0 

[32]  EXTIME1<-C[UIO]x(14^TSLICE) 

[33]  C+-11C  O  C«-“HC 

[34]  D«-0 

[35]  TIMEf^ 

[36]  0 

[37]  n  Average  with  the  value  of  INTV 

[38]  fl 

[39]  INTV«-60 

[40]  If-DIO 

[41]  LIMIT*-(pC]-(DIO+l  +  INTV) 

[42]  n 

[43]  R  Calculate  extra  time  due  to  averaging 

[44]  R 

[45]  EXTIME2«-(  +  /INTVtC)+TSLICE 

[46]  R 

[47]  R  Calculate  the  time  index  which  is  not  equally  placed 

[48]  R 

[49]  LOOP4: 

[50]  TMP«-+/INTVtC 

[51]  D<-D,TMP+INTV 

[52]  TIME«-TIME,C[DIO] 

[53]  C«-11C 

[54]  I«-I  +  l 
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[55]  -♦(ISLIMIT)/L00P4 

[56]  TIME«-(  C  +  \TIME)-5-TSLICE)+EXTIMEl+(EXTIME2-^2) 

[57]  ft 

[58]  ft  Calculate  averaged  frequency  for  the  time  indicies 

[59]  ft 

[60]  Y<-(l-!-2xD)xTSLICE 
V 
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E .  SAMPLE 


[1] 

121 

131 

C4] 

[5] 

[6] 

[7] 

[8] 
C9] 
[10] 
[11] 
[12] 


V  y<-SAMPLE  X;OIO;Y;SIZ;LIMIT;I;X 
o  Routine  to  sample  for  plotting  (number  of  samples 
DIO«-0 

SIZ*-L(  (pX)^600) 

LIMIT<-(pX)+DIO-l 

I«-DI0 

LABELl: 

y«-y,x[ixsiz] 

I4-I  +  1 

-♦(  (IXSIZ)SLIMIT) /LABELl 


V 


60 
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